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The laws of thermodynamics provide a clear concept of the temperature for an equilibrium system 
in the continuum limit. Meanwhile, the equipartition theorem allows one to make a connection 
between the ensemble average of the kinetic energy and the uniform temperature. When a system 
or its environment is far from equilibrium, however, such an association does not necessarily apply. 
In small systems, the regression hypothesis may not even apply. Herein, we show that in small 
nonequilibrium systems, the regression hypothesis still holds though with a generalized definition of 
the temperature. The latter must now be defined for each such manifestation. 
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I. INTRODUCTION 



I A lay person defines temperature as a measure of the heat contained in a body. This is not in conflict with a more 
. rigorous thermodynamic definition based on the equipartition theorem when the body is in the continuum limit with 
respect to size and measurement time. But what is the temperature when these assumptions no longer hold? Is it 
useful to even try to describe small and/or nonequilibrium subsystems using an instantaneous temperature? 

Experimentally, at least, temperature can be defined even for very fast and small systems using heat balances. 
For example, the temperature of a medium can be changed experimentally using a system excitation induced by a 
■ laser pulse or using a microwave field to heat the bath. The temperature in the local environment can be raised 



quickly by several to tens of degrees Centigrade in a few microseconds or lessi This elevated temperature can be 
sustained for milliseconds and subsequently decays back to the global environmental temperature within a time-scale 
of tens of milliseconds Such T-jump experiments provide the time resolution needed for studying the kinetics of 
many processes — as in, e.g., protein- and peptide-folding,~'^>2. and helix-coil transitions in peptides However, 
(~| T-jumps can significantly influence the structural and dynamical properties of the systems under investigation. For 
O ', example, solutions of colloidal and microgel particles change their size dramatically due to changes in temperature 
, — and also pHi^'^ i^^d^i-'^^d^i^'^d^ Indeed, such particles have been swelled by as much as an order of magnitude in just 

I a few milliseconds i^^iii Such environmental changes evidently affect the local and global structure and may enable 
T-H • coupling to additional heat sinks. At the very least, they require additional care in performing the heat balance to 
^ ■ determine the temperature during the T-jump. 

, In order to assign a rigorous definition of the instantaneous temperature of a small (open or closed) system under 
nonequilibrium conditions, it is useful to partition it into a subsystem — specifying the properties of the material — 
and a bath — consisting of everything else. The latter may in turn also be open or closed depending on whether 
or not it is effectively interacting with an additional environment. When this bath is stationary, such as can be 
1^ found in simple liquids or colloidal suspensions at equilibrium, the motion of the subsystem — as represented by an 
n-dimensional position coordinate q — can be accurately described by the generalized Langevin equation (GLE)fti, 



m = - /* 7th(t - mt')dt' + eth(t) . (1.1) 







5—1 ' Here V{q) is the potential of mean force (PMF), 7th(^ — t') is the friction kernel representing the response of the 
I solvent, ^th{t) is the random force due to the medium, and they obey the fluctuation-dissipation relation (FDR), 

(fth(i)6h(i')> =fcBT7th(i-i') • (1.2) 

When the bath is nonstationary by way of exhibiting temporal (and spatial) changes in the ambient bath, Eq. p.l|) 
has been modified to the so-called irreversible GLE (iGLE) !^^'^° 

m = - 1^ lit, t')q{t')dt' + m , (1-3) 



7(t,i') =ffW7th(i-0ff(0 



(1.4) 
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The function g(t) modulating the friction kernel is completely determined by an irreversible process due to processes 
not otherwise included in the subsystem or bath. The stochastic force £^{t) = 5'(i)Cth(i) in the iGLE is modulated by 
the same function g{t). The resulting nonstationary FDR is 

mm) ^ kBT^M ■ (1.5) 

The system-bath dynamics characterized by the GLE and the iGLE have often been generalized to allow for 
the possibility of time-dependent or driven variations in the system temperature. In one set of approaches, the 
nonequilibrium behavior of the bath is introduced by way of an external stochastic force. This additional noise can 
drive the reaction coordinate q directly^^ or by modulation of the local bath modesi ^^'^'^i^'* Although such a modulation 
may be caused by many different physically relevant mechanisms, the recurring requirement is that the external noise 
has no connection to the memory kernel j{t,t') in the FDR. The latter can arise, for example, when the stochastic 
noise acting on the reacting s ubsy stem and the local bath is independent. This requirement necessarily leads to a 
shift in the temperature. Refs. [2lll23l23ll23 have explored this issue within the Kramers theory, and the steady state 
of the nonequilibrium open subsystem was observed to lead to a new effective temperature. 

However, if the system is truly in equilibrium with the overall bath, then the heat transfer leading to this temperature 
differential is presumably dissipated by bath modes not initially included in the local bath. A repartitioning of the 
environment into two sets — a nonequilibrium open reservoir and a global thermally equilibrated bath — have been 
suggested .^^'^^ External perturbations in the initial state of the former are relaxed to equilibrium through the coupling 
to the latter. (The FDR between the random forces and the memory kernel appears naturally in this case.) This 
coupling is stipulated by a specific Hamiltonian model of the bath, though primarily their work has specified the bath 
using the usual choice of a set of harmonic modes coupled by bilinear interactions. Each such mode is then assumed 
to approach its equilibrium value exponentially with a specified dissipation rate. Alternatively, the heat from the 
nonequilibrium local bath may be balanced directly in the iGLE through the ad hoc introduction of an effective heat 
sink which has indeed been seen to provide self-consistent constant-temperature dynamics i^i^Sii^S 

Regardless, the requirement that the overall bath dissipate the system so quickly that constant temperature is 
always maintained is too severe. For example, it does not describe systems in which energy transport fluctuations to 
nearby local bath modes may be important. Nor does it account for (driven) temperature-ramped chemical reactions 
as has earlier been treated using a simple generalization of the iGLE.'^*' That approach is restricted to the case of 
slow temperature variations as will be recapitulated in Sec. IIIIl Thus the question remains as to whether a subsystem 
within a nonequilibrium environment can be characterized using an observable temperature, and if so, how to calculate 
or measure it. 

In the current work, we investigate a generalized chemical process in which the nonequilibrium behavior of the 
chosen subsystem is influenced by a change in temperature of the environment that is in turn influenced by the 
response of the subsystem. In such a case, the environment itself can also persist in a nonequilibrium state. The 
model, therefore, consists of a hierarchy of three domains: a small-dimensional system, a nonequilibrium open "local" 
bath (see Fig. [1]) strongly coupled to the system, and a "global" bath that is coupled to the local bath and at most 
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FIG. 1: A hierarchical representation of a (sub)system under investigation which is immersed within a local nonequilibrium 
environment that is, in turn, interacting with a global equilibrium one. (Color online.) 

weakly coupled to the system. Without loss of generality, the linear response of the local baths can be represented 
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using an auxiliary model of harmonic oscillators coupled bilinearly to the system. The interaction between the local 
and global environments induces a time dependence in the parameters characterizing the local bath — viz., the effective 
masses of the oscillators, their frequencies and the coupling between the local bath and the system. Assuming that 
the effective temperature of each bath mode is observable, we investigate the connections between these variables and 
the system. The time dependence of the bath parameters can be specified according to the behavior of any particular 
physical problem, and so this theory can be applied quite generally as discussed above. 

One important case for these problems (indeed it lies at the core of the discussion that follows) is the possibility 
that the local bath actually consists of several independent baths (or reservoirs). This separation can be viewed 
as a mathematical trick that enables one to demonstrate the possibility of invariances regardless of the specific 
representation. But it is also physically relevant. For example, two- and three- temperature models have been used to 
describe ultra-short laser pulse desorption experiments i'^^'^^ Therein, the distinct reservoirs represent the adsorbate 
layer, the substrate phonons and the substrate electrons, respectively, and the times scales of each can be quite 
different?^ However, these methods have typically been applied in an ad hoc fashion, and the present work is an effort 
towards making this theoretical structure clear in a general sense for cases such as this. 

The GLE can be derived from a Zwanzig-type Hamiltonian as a projection of the simplest many-dimensional 
mechanical systemjiS Within this formalism, the bath is represented as a set of harmonic oscillators coupled with a 
tagged particle by bilinear interactions. In Ref . [20l. this formalism has been extended to take into account nonstationary 
effects by introducing time-defendant coupling coefficients and a nonlocal memory correction term (sec Section . 
The projection of this system has been shown to be the iGLE with the corresponding FDR.'^^ Although the form 
of the iGLE had been known earlier from the thermodynamic considerationsji^ the value of this approach is that it 
allows one to obtain practical stochastic equations at a low cost. For example, the local limit of the iGLE has been 
seen to surmise the probe dynamics in a nonstationary colloidal suspension. '^^ 

In Sec. ini we take advantage of the iGLE and its associated Hamiltonian to derive the stochastic equation for a 
tagged particle moving in a complex nonequilibrium environment consisting of the reservoirs with changing temper- 
atures. We show that this approach can serve as a common basis to unify various methods developed separately 
for specific purposes in Sec. IIIII Moreover, a numerical gedenken experiment looking at the dynamics of a particle 
coupled to three distinct baths (that are not mutually coupled) at different temperatures illustrates the accuracy of 
the nontrivial nonequilibrium effective temperature equation derived in Sec. [Til 

II. LANGEVIN DYNAMICS IN A NONEQUILIBRIUM ENVIRONMENT 

The objective of this work is the determination of the effective temperatures of the "system" shown in Fig. [2l 
The target system is directly coupled to a large number of local reservoirs. These reservoirs are, in turn, coupled 
to a larger scale environment — to wit, the global environment — that alters its properties and interaction with the 
system over much longer time scales than the reservoir relaxation times. In principle, a detailed Hamiltonian can be 
written that would include the degrees of freedom for the system, the local environments and the global environment 
illustrated in Fig. [1] Given such a specification, a projection of the Hamiltonian to the system alone, and a study of 
its correlation functions would reveal its effective temperature and related properties. However, we found earlier that 
the projection onto a system coupled to a single time-dependent local reservoir can be described by the iGLE and 
an associated time-dependent Hamiltonian.^*'''^'' In what follows, we thus take a simpler approach to the solution of 
our objective by generalizing the iGLE and its associated time-dependent Hamiltonian for the case that the system 
can be coupled to many distinct local reservoirs. While much of the discussion takes advantage of the form of the 
auxiliary Hamiltonian system to derive various correlation functions, all of the latter involve only variables of the 
reduced dimensional system. As such, these results are independent of the particular specification of the auxiliary 
Hamiltonian. 



A. Multiple-reservoir model with fixed temperature 

The Brownian motion of a tagged — or chosen — particle immersed in a single equilibrium thermal bath is well 
understood. More generally, however, the local environment of the particle can be further divided into a set of 
different reservoirs — viz. distinct local baths — with separately identifiable characteristic properties. In this section, 
all of these reservoirs are assumed to be at the same temperature and therefore independently satisfy the same 
Boltzmann distribution over their coordinates, (c./.. Fig. [3]). A more general case in which these reservoirs may follow 
different temperatures — and consequently, as a whole, the environment does not follow a simple sing-temperature 
Boltzmann distribution — is addressed in the next section. 
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FIG. 2: A schematic representation of the generalized model described in this work. The subsystem interacts with multiple 
bath reservoirs by means of time-dependent couplings while the temperature in each bath is also subject to change. (Color 
online.) 



W) 




System 



Ut) 




FIG. 3: A schematic representation of the model investigated in Refs. [l9l. \2ui and ISa. The subsystem interacts with multiple 
bath reservoirs by means of time-dependent couplings while the temperature remains constant. (Color online.) 

In order to describe the dynamics of a particle solvated by many distinct reservoirs, one can generalize slightly the 
Zwanzig-type Hamiltonianj^i^ 



H 



^ + Viq) + SVi{q,t) + SV2{q{-),t) 



(fc) 
b ' 



(2.1) 
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where each reservoir, fc, is represented by a distinct set of harmonic modes. The modes in each reservoir interact 
directly only with each other and the tagged particle, and indirectly to modes in other reservoirs through their mutual 
couphng to the tagged particle. In the Hamiltonian of Eq. p.ip . the sum of the first and second terms, + V{q), 
is the bare energy of the tagged particle. The last term is a sum of Hamiltonians, 



Hi"' = 



Pi 2;, 



(2.2) 



(k) (k) 

each of which represents the energy of the fc-th bath reservoir with coordinates and momenta pl . The fifth term 
provides the bilinear coupling between the bath modes and the particle, in which and denote the summation 
over the bath reservoirs and the bath modes (i.e., bath frequencies). The third term, 

1 / r^^') \ ' 

^V^i(9,t) = ^EE3feWHfe) (2-3) 

hi / 

is the renormalization of the potential which eliminates the time-dependent spectral shift, and the fourth term, 

5V2{q{-),t) = \ f dt'a{t,t')[q(t')-q{t)f-lq{tf f dt'a{t,t') , (2.4) 
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provides the time-dependent correction to the memory effects.— 1^ The nonstationary memory correction in Eq. (|2.4p 
was defined earlier as 

a{t, t') = 9k{t)gk{t'hith)k{t - t') , (2.5) 

k 

where 

/ \ ^ 

7(th)fc(t -i') = Y.\jr)] cosa^f'V - t') (2.6) 



is the friction kernel for the fc-th reservoir at thermal equilibrium. 

The nonlocality in the 5V2{q{-), t) term has been addressed earlier^^iJS and is necessary because of the time depen- 
dence in the coupling terms. (If said time dependence were zero then this term would also be zero.) The former, in 
principle, also arises in the case of space-dependent frictio n'^^i^^i'^^i^^i"^° — viz. space-dependent coupling — where the 
factors gk in the fifth term of Eq. I|2.ip depend only on the coordinate q. Although the transient coupling terms have 
recently been noted in said case^i^^ they have been completely discarded in the corresponding Hamiltonian equations 
of motion by way of heuristic arguments concerned with the possibility that they appear to violate causality. In the 
present case, however, where several coupling terms depend explicitly on time, the nonstationary memory correction 
SV2{q{-),t) — requiring the integration over the particle's path g(-) over all time — cannot be discarded. Indeed, this 
correction is necessary so as to ensure that the projection onto the iGLE does not itself contain additional nontrivial 
terms as detailed in Appendix [Bl that would violate causality. As detailed earlier the component of the nonstationary 
memory correction that violates causality does not present any problems because it is smaller than other terms that 
have also been neglected in this formalism. 

In Ref. [2^, the iGLE — Eq. (|1.3p — (with a family of nonstationary friction kernels) was shown to be the projection 
of the equations of motion derived from Hamiltonian (j2.1[) when there is only one reservoir. In the current case 
of multiple reservoirs, the friction kernel in the iGLE can be represented as a sum of equilibrium friction kernels 
modulated by distinct functions gkit) defining the coupling strength to the corresponding bath reservoirs, 

lit, t') = 9kit)gkit'hitb)kit - t') . (2.7) 



The functions gk (t) are specified by the time-dependence of the coupling up to a trivial arbitrariness in their initial 
values that is renormalizable. In cases when the reservoir is initially uncoupled to the particle, then said values are 
necessarily specified and equal to zero, however. 
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The FDR (jl.Sp for the projection of the iGLE, Eq. (|1.3p . with the multi-reservoir friction, Eq. (|2.7p . results from 
the Hamiltonian of Eq. (|2.ip , if the initial distribution is of the form, 



P oc cxp 



k,i ^ 



ujI 'xl '{0) 



UJ, 



(k) 



9(0) 



and the conditions 



(2.8) 



(2.9a) 



„(0 



5;(0) 



,(0 



9(0) 



= kBTSijSki , 

/ i'')/n\ ( (0 (0/n^ 4 ^ff'(O) 



9(0) 







(2.9b) 



(2.9c) 



are satisfied. 

The iGLE with multiple and distinct reservoirs has thus been specified, and allows for temporal changes in the 
response of the environment albeit while all such reservoirs maintain the same constant temperature. 

B. Multiple-reservoir model with time-dependent response 

1. Squeezing bath modes 

In principle, the effective temperature of each of the reservoirs (solvating the tagged particle) need not be the same 
as they are not coupled to each other directly. This situation is illustrated in Fig. [2] in which the temperature of 
each reservoir is different, perhaps as a consequence of some unspecified external or internal processes. It can be 
implemented within the Hamiltonian representation of a particle interacting with multiple reservoirs by supposing 
that the bath modes are not static — that is, their effective masses and coordinate scales are allowed to vary in time. 

The time-dependent functions, /ife(i) and /i/t(t), are now introduced so as to represent the time-dependent scaling 
in the force constant and mass of the fc'^ bath mode, respectively. The Hamiltonian of the i-th oscillatory bath mode 
belonging to the /c-th reservoir can thereby be written as 



2/4(0 



(2.10) 



The terms, hk{t) and Hklt), combine to determine the time dependence of the mode frequencies in the fc-th reservoir. 
Their initial values can be chosen as /ifc(O) = /ifc(O) = 1 without loss of generality. Meanwhile, the scaling of the force 
constant effectively narrows (or widens) the configuration space available to the oscillator at a given temperature. 
Such a trick of "squeezing" bath oscillators leads to changes in the total energy of the reservoir, and thereby shifts 
the reservoir temperature. The existence of such temperatures relies on the fact that the relaxational redistribution 
of the collective energy within the modes is sufficiently fast that bath reequilibrates on the time scales of the particle 
motion. In other words, the process of bath mode squeezing is assumed to be adiabatic and to vary on time scales 
longer than the periods of oscillation in the nontrivial modes of a given reservoir. 
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The time-dependent Hamiltonian (|2.ip can now be modified by introducing the scahng functions hk(t) and fJ.k{t)'. 

2 

H - ^ + Viq) + SViiq,t)+SV2{qi-),t) 

k,i k,i 

(2.11) 

Here the coupUng coefficients 5^ (t) are introduced so as to simpfify the structure of the Hamihonian. These coefficients 
are unequal to the earher coefficients gk{t)- Nevertheless they can be connected through the expression, 



9'kit)^gk{t).H§-, (2.12) 

allowing comparison to the earlier results obtained using the original iGLE i^^'^°i'^°'^^ In Eq. (|2.1ip the renormalization 
potential, SVi , is given by Eq. (j2.3p with gk (t) replaced by g'j^ (t) . The nonstationary memory correction, 6V2 , is defined 
by Eq. with 

«(^'^') = E ^.9feW3fc(^')7(th)fc [Tk{t) - n{t')] , (2.13a) 
where gkit) has not been replaced by g'^it), 

Tkit)= f hk{s)|^ik{s)ds (2.13b) 
Jo 

and 7(th)fc(') is defined in Eq. (|2.6p . 

Note that the coupling coefficients gk{t) can also depend on the temperature. Although in the framework of this 
approach, there is no recipe for connecting the temperature-driving functions /ifc(i) and hk{t) to these coefficients, 
such a connection can in principle be obtained from a deeper analysis based on mode-coupling theory, direct extraction 
of these functions from simulations (Ref. i35i: see also Sec. IIIIC] below), or experiments. 

As was done in the earlier section with respect to the behavior of gk{t), the functions fJ.k{t) and hk{t) must be 
required to be effectively stationary on the correlation time scales of the fastest bath modesJ^i2£ This so-called 
slow-bath-mode "squeezing" assumption is necessary because otherwise there would always exist a bath mode whose 
response to fluctuations would be so chaotic that it wouldn't ever relax to quasi-equilibrium limits. Evidently, for low 
spectral frequencies, this imposes some restrictions on the rates of change of fikit) and hk{t) and, thus, implies a related 
assumption that low frequency bath modes contribute negligibly to the particle's motion. These assumptions are valid, 
for example, in the over-damped case when the spectral density, J{oj), of the bath Hamiltonian is proportional to 
with s > 1 at w — > 0, and hence the small bath frequencies are not coupled to the system. 



2. The effective temperature of each reservoir 

The strategy for solving the dynamics of Eq. (|2.10p . taking into account its hierarchical structure, focuses first on 
the solution of the local dynamics of the bath modes within a given fc'^ reservoir as if it does not interact with the 
rest of the modes. In the next section, we will use these states as a reference for the approximate solution of the 
global dynamics. The accuracy of this approach relies on the assumption that the coupling between the reservoirs 
and the tagged particle is weak, as has indeed been assumed from the outset. For simplicity, in this subsection the 
upper indices (fc) will often be omitted as all of the formulas herein are necessarily being solved within a given fc*** 
reservoir. 

At any given time, the energy manifold of the i*^ bath mode in the fc**^ reservoir is simply a function of the phase 
space variables, {xi,pi), determined by the instantaneous parameters, i.e., 

m)^^^ + d,,^tr,^^. (2.14) 

If Hkit) and hk{t) change slowly relative to the oscillator frequency, then the symplectic area of this manifold — viz., 
the action / — will be an adiabatic invariant. On the elliptical energy manifold defined by Eq. (|2.14p . this invariant 
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is given by 



hk{t)u 



hk(0)uji 



(2.15) 



where the last equahty foUows from the initial choice, hk{0) 
useful connection, 



ES) = E,{0) 



fJ-kiO) = 1- The time-invariance of / thus provides a 



hk{t) 



(2.16) 



between the energies of the i^^ bath mode at various times. 

Given that each bath mode is initially equilibrated at its corresponding reservoir temperature Tfc(O) — i.e., the 
initial energies £"1(0) of the reservoir are Boltzmann distributed, — the connection in Eq. (|2.16|) readily implies that 
the energies Ei{t) at any given time are also Boltzmann distributed at the temperature, 



Tk{t) = r,(o) 



hk{t) 
M/c(i) 



(2.17) 



Thus, the ratio hk{t)/iik{t) establishes the time dependence in the temperature of the fc*'' bath reservoir. Although 
the interaction between the reservoir and the tagged particle could, in principle, alter this structure by way of heat 
transfers, in practice it does not do so. The bath reservoir is assumed to be sufficiently large that such transfers are 
small in comparison with the total energy of the reservoir. (This assumption is not very severe because it must be 
satisfied in order for the bath to have a well-defined quasi-equilibrium temperature!) The functions hk{t) and fik{t) 
are auxiliary and are not known a priori, but the time-dependent temperature entering in Eq. (|2.17|) . together with 
the coupling strength coeffrcients gk{t), can either be found from simulations or estimated from experimental data. 
For example, they have been obtained in the case of a diffusing particle in a swelling colloidal suspension in Ref. [ssl . 

The initial temperatures of the various bath reservoirs need not be the same because some modes may be con- 
nected to distinct heat sources and sinks leading to steady state energy transport through the system between the 
reservoirs. Indeed, the nonequilibrium dynamics of a particle connected to reservoirs with different time-independent 
temperatures has been explored by Kurchan and coworkers j ^^'^*^ But if the system is initially in equilibrium, then 
in the absence of such energy flows, the initial temperatures of all the reservoirs will be the same. More generally, 
however, the initial probability distributions — cf. Eq. (|2.8p — for a given reservoir can be rewritten as. 



Pk oc exp < - 



2fcBrfc(0) 



-f^.f(0)-^Jff^^.(0) 



(2.18) 



for the fc*^ reservoir temperature, Tk{Q). The correlations in Eqs. ()2.9a|) - ()2.9cp are still satisfied if one replaces the 
temperature, T with that of the given reservoir, Tfc(O). Such an initial equilibration seems to be artificial when the 
temperatures Tk change with time due to the interactions of reservoirs with the tagged subsystem and the global 
bath. Note, however, that the initial time, t = 0, can be formally shifted to the far past. The system then has no 
memory of the initial conditions and, thus, the resulting iGLE remains unaffected. 
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3. Equations of motion 



The equations of motion for the many-reservoir case foUows directly from the Hamiltonian in Eq. (|2.1ip : 



{k) 



.(k) 

q 

Pi 



Pg 



k.i 



d 



(k) 



(2.19a) 

(2.19b) 
(2.19c) 
(2.19d) 



dqit) 



SV2{q{-),t) 



where the indices, k, referring to a given reservoir have been exphcitly retained to help the reader. However, they 
are somewhat cumbersome, and so in this section it is convenient to replace explicit reference to the superscripted 
reservoir index. Instead the collective index i(= {i, k}) will be used for this reference, as in e.g. coi = w^^. 

As discussed in Sec. Ill A| the presence of the nonstationary memory correction d{6V2)/dq{t) in Eq. ()2.19|) introduces 
the possibility of violations of the causality principal. Indeed, the variational principle applied to the calculation of 
this term gives an apparent contribution to the force from future trajectories. However the error in this contribution 
is on the order of the neglected terms in the perturbative treatment in the projection, and it was seen earlier'^^ that 
all of these small errors can be safely neglected. But the nonstationary memory correction d{6V2) / dq{t) can not be 
entirely neglected because it contains a nontrivial correction which ensures that the projection to the iGLE docs not 
violate causality up to the second order of the perturbation expansion (c.f.. Appendix IB|) . 

From Eqs. (|2.19ap and (|2.19bp one can derive 



ii + i^kit)ii + il?(i)xi 



„ 9'k(t)hk{t) ^ 



where 



iyk(t) 



2Afc(^) 
hkjt) ^ 

^J.k{t) 



(2.20) 

(2.21a) 
(2.21b) 



The assumption of a slow "squeezing" bath modes, discussed before, must now be restated in terms of both time 
dependent parameters, hk{t) and fj.k{t). After a little algebra, it can be shown to reduce to the requirement 



(2.22) 



which must be satisfied for all modes i in a fc"^ reservoir. Under this assumption, one obtains the following result [as 
is shown in Appendix [X)) ] : 



Xi = Ui{t)xi{0) + Ui{t)xi{0) 



dt' 



CiWi(i) 



(^i(i')Mfc(i'))' 

dt' 



{u,{t')^iu{t')Y Jo 



g'f,(s)hk{s)ui{s)q{s)ds , 



where the bath mode oscillations are described by 



Ui{t) 



cos LO\Tk{t) 

^tik{t)hk{t) 



(2.23) 



(2.24) 
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and the effective time can be written as 



Tkit) 



1 



n;{t')dt' 



hkjt') 
f^k{t') 



dt' 



TkiO) 



dt' 



(2.25) 



Perhaps not surprisingly, the effective time changes "faster" when the temperature increases. A similar transformation 
of time has been found earlier in simplifying overdamped Langevin equations driven by time-dependent temperature 
baths4^^^ 

The projected equation of motion for the tagged particle reduces to the iGLE form after the appropriate averaging 
over the distributions derived in the previous subsection. (The details of this derivation can be found in Appendix [Bl) 
The nonstationary friction has a form similar to that of Eq. (jl.3p with 



l{t,t') 



E ^fk{t)gkit'h^ty.)kirkit) nit')) 



for t > t', where the weighted friction of the k mode is defined as 



lkit,t 



,^ _ ^Tu{t)Tu{t') 
Tk{0) 



9k{t)gk{t'h{th)k{Tk{t) - Tk(t')) 



(2.26a) 
(2.26b) 

(2.26c) 



and is symmetric in the two times, t and t'. (Note that 7(th)fe(') is an even function of its argument as can be seen in 
Eq. (|2.6p .) Recalling the definition for the bath reservoir temperature and friction kernel in Eqs. (I2.17|) . and Eq. (|2.6|) . 
respectively, the nonstationary stochastic force is found to be 



9k[t)y — 



k ^'=(0) 



Pi(0) sinwiTfe(t) 



^.^.(0)-£iM^g(0) )cosc^irfc(i) 



Imposing the initial condition in Eq. (j2.18p . the random force correlation function becomes 

Tkit)Tkit') 



mm) 



k 



Tfe(O) 



-9k{t)gk{t'h{th)k(Tk{t) - Tfe(t')) 



= a:b E VnWW)ikit, t') . 



(2.27) 

(2.28a) 
(2.28b) 



Equations (|2.26p and (|2.28p can be used to recast the extended form of the FDR, 

m)m)) = kBT,n{t.t'h{t,t') , 

where the effective temperature is 



Tcff(i,t') 



Y.^,^Tk{t)n{t')-ik{t.t') 



(2.29a) 



(2.29b) 



for t > t' . Eqs. ()2.29p comprise the central result of this work and generalize the FDR by including the simultaneous 
interactions of the tagged particle with different nonstationary thermal reservoirs. In particular Eq. (j2.29b[) provides 
a clear prescription of the effective temperature of the chosen particle that emerges as a consequence of averaging 
the microscopic dynamics of many reservoirs. This is a nontrivial result because it dances between microscopic and 
macroscopic quantities. Moreover, as written, the effective temperature appears to be asymmetric between the two 
times, t and t', and hence may cause alarm. However, only the product Toff and 7 in Eq. (j2.29ap must remain 
symmetric. As shown in the discussion that follows, this product is indeed symmetric with respect to the two times 
in its argument, and T^e reduces to forms that have earlier been derived or assumed for a number of limiting cases. 

One possible concern with the underlying equations that have led to the Eq. (|2.29p lies in the specification of the 
initial condition for the system and its environment, given that there must evidently be some environmental memory 
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of the previous motion by way of the values of q(t) at earher times. In past and present work, we have typically 
employed one of three boundary conditions: (i) The system is at equilibrium at the initial time, 0, just before the 
nonequilibrium disturbance is effected. This requires us to run dynamics starting at a time well prior to t = allowing 
the system to equilibrate and fully prepare its environmental memory. (This evidently requires sampling over many 
trajectories.) (ii) Different reservoirs att = are at distinct equilibrium within themselves, but are not at equilibrium 
with each other. In this case, we prepare each reservoir separately using the procedure described for case i. (Hi) The 
particle is suddenly placed inside the reservoirs, and hence there is no prior memory due to motion before t = 0. In 
this section, we have applied the boundary conditions from case ii, in order to simplify the expressions from correction 
terms that would have been trivial anyway. Indeed, the use of case Hi for the boundary conditions in the numerical 
simulations shown below did not affect the results. 

It may appear that the terms in the central result of Eq. (|2.29p are contingent on the creation of an associated 
Hamiltonian of the form of Eq. (j2.10p used in its analytical derivation. However, as in other such reconstructions using 
the latter auxiliary infinite-dimensional representation, the projected form of the central result contains no object with 
explicit reference to any unprojected dynamics. As such, it should be capable of describing any dynamics which is 
captured by the continuum representation of the iGLE. More specifically, the temperatures Tk can be calculated or 
measured directly from each of the independent baths. The coupling terms gk and memory times Tk can be calculated 
or measured directly by obtaining the correlation functions of the forces exerted on the tagged particle by each bath, 
respectively, upon removal of the thermal component. Such a procedure has been used in Ref. 35 in the limit in which 
the chosen particle is coupled to only one bath. Several additional examples are provided below in which the particle 
is coupled to multiple baths. In each of these cases, the associated functions, gk{t), Tk{t), and Tk{t), necessary for the 
application of Eq. (|2.29p are constructed. It is notable, that these functions can be obtained without recourse to the 
auxiliary infinite-dimensional form of Eq. (j2.10p . 

III. RESULTS AND DISCUSSION OF LIMITING CASES 

In order to better understand the effective temperature expression derived above, it is helpful to look at a number 
of limiting cases. In Sec. IIII Al the structure of the correlation functions when all the reservoirs are held at constant 
but distinct temperatures is shown to agree with the earlier perspective of Kurchan and coworkers4i In Sec. IIIIBi 
the generalized Langevin equations used earlier to describe aging systems are recovered for when the reservoirs are 
all held at a homogeneous, but time-dependent, temperature. This latter limit is also equivalent to that suggested 
earlier in Ref. |35| in a case studying driven nonequilibrium colloidal suspensions. In Sec. IIII C 21 the two-temperature 
model of Refs. [25| and 26 is also recovered in the limit that the system is coupled to two local reservoirs. A set 
of gedenken experiments wherein the bath is coupled to three local reservoirs is shown in Sec. IIII C 31 to lead to a 
nontrivial effective temperature in agreement with the central results of this work. In Sec. IIII C 41 the local quasi 
equilibrium limit of a non-Markovian Fokker-Planck equation explored earlier by Rubi and coworkers^ is shown to 
lead to correlation functions equivalent to those in the corresponding limit of this work. Thus the central result for the 
effective temperature of the system driven by nonequilibrium local environments is shown to capture all the previous 
limits as well as extend the formalism to heretofore unknown cases. 

Before proceeding, it is perhaps also useful to recapitulate the separation of time scales in the hierarchy illustrated 
by Fig. [T] that has been assumed in constructing the central result in Eq. (|2.29p . The relaxation times for responses 
to large perturbations of the system Tg, the local reservoirs t^, and the global environment Tg must satisfy the simple 
ordering, Tg < < Tg. Meanwhile, the system is so small compared to the local reservoir that any perturbation 
at times on the order of Tg leads to no effective change of the local reservoir behavior. Any sustained change in the 

,th 

system behavior over times longer than r^,' will affect and change the behavior of the k reservoir, but it will simply 
lead to a change of the quasi-static response of said reservoir to the system at the short times scales near Tg. An 
analogous description is true for the time scales between the local reservoirs and the global bath. Because of these 
inequalities, any time-dependence in Tk and gk must be slow in comparison with the motion of the system. However, 
as we have seen in earlier work, and shown here in the numerical work, the time scales do not have to be nearly so 
disparate as the corrections are small. 

A. Environments with Constant Temperature Inhomogeneities 

The structure of Eq. (|2.28p and the associated stochastic equations derived from them allows one to describe the 
dynamics of a chosen particle connected to many reservoirs where each can be at a different temperature. If these 
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temperatures are constant in time, i.e., Tk{t) = Tk{t') = Tk, then from Eqs. (|2.25p . (|2.26ap and (|2.28p one obtains 

7(f,t') = Y.9k{t)9kit'hitb)kit~t') , (3.1) 

k 

mat')) = 5]M5fe(i)5fe(<')7(th)fe(t-t') ■ (3.2) 

k 

An analogous case has been discussed in Ref. in the stationary Hmit when (t) = 1 (stable non- uniform environ- 
ment). Thus, Eqs. p.ip and (13. 2p extend the previous treatment so as to include nonequilibrium baths. 

The implementation of the iGLE theoretical framework discussed above follows readily for any given nonequilib- 
rium system interacting with many nonequilibrium baths. The key step is the identification of which baths — viz., 
temperature heat sinks — are in contact with the chosen particle as a function of time. Nonzero and zero values of 
gk{t) at any given time correspond to whether said reservoir is connected and disconnected, respectively, from the 
particle at the given time, t. (An analogous procedure was used in Ref. 49, for deriving a variant of the fluctuation 
theorem.) This will help to switch the chosen subsystem between different heat reservoirs with different temperatures. 
If the reservoir temperature does not change significantly during a period when an appropriate coupling term differs 
from zero, one obtains Eqs. (|3.ip and (|3.2p . 



B. Environments with Time-Dependent Homogeneous Temperature 

If the temperature changes uniformly for all the bath modes, Tfc(t) = T{t), then Eq. (I2.28P simplifies to 

mm) = kBT{t')^{t,t') . (3.3) 

Such generalized FDRs (GFDRs) have been used in the description of the nonequilibrium steady-state dynamics of 
glassy systems, for examplc j^?'' ^- The relations between autocorrelation and response functions, like that of Eq. (|3.3p . 
are frequently considered in the literature as a signature of FDR violations (quasi-FDR) arising from partial equili- 
bration among a subset of degrees of freedom, and the quantity T{t') is treated as an effective temperature. In the 
specific case of aging systemS )^°'^^i^^ the effective temperature depends on both t and t' . The name "quasi-FDR" and 
the use of T^s as an auxiliary parameter are reasonable due to the fact that the actual thermodynamic temperature 
is constant in these studies. 

The results of the previous section, however, allow for the recognition that T(t') is the true ambient temperature 
at the time t' , and the interpretation that GFDR in Eq. (|3.3p expresses the transient behavior of the system. Note, 
that in the iGLE model, the partial relaxation of the baths corresponds to different dependencies Tk{t) [included into 
the GFDR, Eq. (H^H])], and leads to the two-time dependence of T^ff in Eq. (|2.29ap . 

In the specific case when all the coupling coefficients are also equal, 5fe(t) = g(t), we obtain 

lit, t') = |||5(t)g(t')7th(r(t) - T{t')) , (3.4) 
for t > t' , where 7th (t ^ '''') = X^fc 7(th)fe(''" ~ ''"') is the friction kernel of the whole system at equilibrium, and 

m 

T{0) 



T(t) 

m - T^.amti^irit)) , (3.5a) 



with 

(^th(T)6h(T')) = A:BT(0)7th(T - t') . (3.5b) 

It is instructive to compare these results with those obtained in Ref. [s^ from phenomenological considerations. 
Namely, the friction kernel used in that work, 

7(t,i')-5(t)5(07th(i~t'), (3.6) 

coincides with Eq. (|1.4p and may include effects due to the temperature change implicitly through the coupling 
coefficient g{t) only. 

This restriction on the memory kernel also induces the main difference between the two approaches: the use of 
the real time t vs. the effective one, T{t). Hence, only the amplitude of the stochastic force changes along with the 
temperature alteration, and gives 



^(i) = W||y5(i)Cth(0. (3.7) 
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The use of the real time instead of the parameterized one imphes that the random force autocorrelation function 
acquires a different form: 

mm) - kBVmW)9it)9{t'hthit - f) = k^vrnWhit^t') . (3.8) 

This FDR varies from Eq. ()3.3|) insignificantly if the temperature changes slowly during the solvent response time, 
when 7th(^ — t') differs from zero. However, the use of the effective time to describe the behavior of the Gaussian 
noise in Eq. p.Sp provides an important contribution that allows the correct projected equations of motion to better 
control the intensity of the stochastic force when the bath is far from equilibrium. This is the main difference between 
the Hamiltonian approach of the current work and that of Ref. i3^ where the random force is governed only by the 
changing amplitude. 



C. Stochastic Dynamics 



The friction kernel of a free heavy particle, diffusing in solution without external potentials (V{q) — 0), has a 
memory-less form^ 

7(th)fc(T) = 27ofc(5(T) . (3.9) 

Recalling the standard functional relation, d{f{x)) — 6{x — a;o)/|/'(a;o)| j where for simplicity we have assumed that 
xq is the only root of the argument, i.e. f{xo) = 0, then Eqs. (|2.26ap and (|2.28p lead to 

7(t, t') = ^ glit) ■ 2^okS{t - t') = 27j{t)S{t - t') , (3.10) 

A; 

{mW)) = fcB ^ n{t)gl{t) ■ 2jokS{t - t') EE fceTeff (t)7(i, t') , (3.11) 
where the effective temperature, Teff (t), of the Brownian particle is set to 

Joff(i) = — — 27T\ ' K-i-^^) 

Ekdiithok 

so as to satisfy the FDR. 

The effective temperature is consistent with that of Refs. andlH for the special case when the noise is memoryless 
and all the reservoirs are stable in time — i.e., when the temperatures Tk are constant and gk = 1. Herein, the resulting 
iGLE turns into the memoryless irreversible Langevin equation (iLE): 

m = -vit)m+m , (3.13a) 

m^it')) = 2kBT,s{t)Tjit)S{t ~ t') . (3.13b) 

The form of this stochastic equation of motion is similar to that of a particle diffusing in a medium with spatial 
temperature gradients4^ 

q = -viq)q + mt) , (3.i4a) 

mmq^t')) = 2kBT{qUq)6it-t') . (3.14b) 

In both cases, the particle quickly equilibrates as it traverses from one local region to another, in conformity with the 
memoryless limit. Such a traversal is reminiscent of space-dependent friction dynamics wherein the stochastic particle 
traverses between distinct dissipating environments. The difference here is that the dissipation can change even if the 
particles stays in the same region because of changes in the dissipating environment. Thus the earlier Darboux-like 
topological construction of the dissipating environments is now extended to be time-dependent. 



1. Uniform temperature Limit 



In the limit that all the reservoirs obey a homogeneous temperature profile, Tk{t) = T{t) for all fc, so that Toff(t) = 
T{t) (see Eq. p.l2p ). the friction coefficient in Eq. (|3.13p becomes 



r;(i) = G2(i)7o, 



(3.15) 
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where G{t) is the effective couphng coefScient defined according to Ref. [H, 

G^{t) = 5Z dkithok/ 70 and 7o = • (^-l^) 

We thus recover the iGLE of Ref. [H in which the stochastic particle is effectively coupled to a single homogeneous 
(but time-dependent) reservoir. 



2. Two-Reservoir Limit 

In Refs. and [1^, the LE for the Brownian particle is derived in the case that the latter is dissipated by two 
reservoirs: a global thermal bath at temperature T, and a local nonequilibrium (time-dependent) bath. The LE 
reads, as usual, 

m = -(r/cq + VncMt) + C(0 > (3-17) 

but with the stochastic force subdivided into equilibrium and nonequilibrium parts, ^{t) = ^cq{t) +S,ncq{t)- The FDR 
for the equilibrium part is (^eq(i)Ceq(i')) = '^VeqkQT6{t — t'). The energy density of the nonequilibrium bath modes 
changes with time and was found in Refs. '25' andll^ to be 

uiuj,t) = / dr(^„eq(Oencq(i + T))e^'^^ (3.18a) 

4?7noq J-oo 



kBT 



"('^-O) — 



(3.18b) 



This takes into account the average dissipation — v in the notation of the earlier articles — of the nonequilibrium 
reservoir modes due to their coupling to the thermal reservoir. 

Within the framework of the present approach, all the modes of the nonequilibrium reservoir have one and the same 
initial temperature, u{uj,0) = kBTo/2. The time-dependent temperature of this reservoir according to Eq. (jS.lSp is 

Tncq(t) = T + e-^'/^To - T) , (3.19) 

and it is the same for all the modes. Taking the inverse Fourier transform of Eq. (|3.18p . one obtains 

(e(i)C(t + t)) = 2kBT,ff{t)rj5iT) , (3.20) 

where 

Teff(t) = '?°q^ + '?"°q^"°qW . (3.2i) 

^cq ~t- ^neq 

This latter result is in agreement with the effective temperature in Eq. (|3.12p if gk{t) = 1. The latter requirement is 
precisely what is needed to ensure that all of the nonequilibrium baths reduce to the single local nonequilibrium bath 
under consideration in this subsection. 

Hence the current construction of the generalized effective temperature of a particle connected to an arbitrary 
number of baths reduces to the appropriate limit when it is connected to two nearly-separable baths. 



3. Three-Reservoir Gedenken Experiments 



To further illustrate the accuracy of Eq. (j3.12|) in providing the effective temperature of a subsystem interacting 
with nonequilibrium reservoirs, we construct a simple model for the diffusion of a heavy spherical particle, with mass 
M, immersed in a three-component solvent gas of light spherical particles. The particles in each of the components 
are assumed to have the same mass m. Each component particle interacts only with other particles of the same 
component through hard-sphere collisions so as to quickly achieve a quasi-equilibrium temperature Tk for the k^^ 
component. All the particles interact with the heavy particle through hard-sphere collisions. A critical assumption is 
that each component does not interact with the other component, and thus each can maintain a distinct temperature 
for very long times. This severe assumption is the origin of our use of the term "gedenken" to describe this model. 
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In principle, a numerical simulation of this system can be accomplished using a many-particle Molecular Dynamics 
simulation in which each solvent component does not in any way interact with each other. However, if the desired 
observables are restricted to correlation functions at times longer than the mean time to collision, then these assump- 
tions are sufficiently severe that the dynamics of the system can instead be performed using a statistical approach 
ala the Enskog theory for gases. Briefly, the algorithm involves the free propagation of the hard sphere particle in 
three-dimensional space for a duration that extends to a randomized time that is consistent with the statistics for 
collisions with each of the particles. The collision takes place with a particle of type k whose 3-dimensional momen- 
tum is consistent with a Boltzmann distribution at temperature T^. After the collision, the heavy mass has a new 
momentum and proceeds to the next collision. 

More precisely, the time tk to the next collision between the heavy particle and a particle from the k-th reservoir 
is calculated as a Poissonian stochastic variable with the average value 

tk = , (3.22) 

where R is the radius of the solute particle, is the radius of particle in the fc-th reservoir, and is the number 
density of the fc-th reservoir. The mean relative velocity Uk{V){= {\V — Vk\)) between the Brownian particle velocity 
V and the particle velocity Vk of the fc-th reservoir is given in Eq. (|3.24p below. In the numerical simulation to be 
discussed, only three reservoirs are included — that is, fc G {1,2, 3} — but the algorithm would still be efficient for an 
arbitrary number of reservoirs as long as there are finite in number. Given the stochastic times, ti, t2, and ts, the 
heavy particle collides with the reservoir particle whose corresponding collision time is the shortest. The time t* and 
label fc* of the closest collision corresponds to this minimal . 

To calculate u^, we again invoke the assumption that each fc*'' reservoir is distributed according to a Boltzmann 
distribution that is unperturbed by its small interaction with the heavy particle.— Thus, the Maxwellian distribution 
of particles of the fc-th reservoir, 

„ s f Tn \ / mv?. \ 



leads to the average relative velocity,- 
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2kBTkJ V V V 2fcBr; 

where V is the velocity of the heavy solute particle and the error function is defined as erf (x) = 2/y/n- cxp{—y^)dy. 

Before the collision, the colliding bath particle is assigned a vector w^. whose components are calculated from 
the corresponding Maxwellian distribution (I3.23p . After the collision, the change of the velocity of the heavy solute 
particle is defined as 



dV 



M 



(vk. ~V+\vk. ~V\-nj 



where n is the unit scattering vector. This vector is distributed isotropically when hard spheres collide. 

Three different gedenken simulations have been performed with the parameters prescribed as follows. The volume 
and the mass of the solute particle are taken to be 64 times larger than that of the reservoir particles which suffices 
to ensure the assumed time scale separation. The specific values for the masses and radii are taken to be M = 
5.181 • 10^^^ g, m = 8.095 • 10~^° g, i? = 4 nm, and = 1 nm. This is consistent with the known density of bulk 
gold — ca. 19.32g/cm'^ — and corresponds to large and small spherical clusters with about 16,000 and 200 gold atoms, 
respectively. Such clusters have been readily constructed in the literaturCf^i^ though of course the separation of the 
light clusters into three distinct reservoirs would be somewhat harder to maintain for extended times. 

The number densities of the reservoirs are taken to be pk = 1.919 • 10~^ nm~^ corresponding to a dilute gas at 
pressures on the order of 10"** Atmosphere. The temperatures of the first and third reservoirs are Ti = 100 K and 
Ta = 500 K, respective. The temperature T2 of the second reservoir is varied in each of the simulations as noted in 
TablelH As the number densities of the three reservoirs are equal, their average temperature, Tavg = (Ti + T2 -|-T3)/3, 
might naively be expected to be the effective temperature of the heavy particle. In order to not bias the results away 
from this naive estimate, the heavy particles are initialized with velocities drawn from a Maxwellian distribution at 
Tavg- We found that good convergence and adequate sampling for the correlation functions of interest to this work 
could been obtained for each case by sampling 10,000 trajectories for up to 4ms. 

The results of the simulations are presented in TablelH where Ty is the relaxation time of the velocity autocorrelation 
function (VACF) of the heavy solute, te is the relaxation time of the solute's kinetic energy, and rj'™(= -^{V^)) is 
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T2, K 




relaxation time, 






temperature, K 








2rE 


rpsim 
Jeff 


TntflCO 

-'off 


Tavg 


100 


0.3513 


0.3483 


0.3621 


307.1 


309.3 


233.3 


300 


0.2984 


0.2978 


0.2992 


346.4 


348.6 


300.0 


500 


0.2704 


0.2703 


0.2701 


422.3 


425.4 


366.7 



TABLE I: The results of a numerical gedenken experiment are summarized in this table. The chosen particle is coupled to 
three distinct baths with temperatures, Ti, T2, and T3, respectively. In all three experiments, Ti{= lOOK) and T3(= 500K) are 
held constant while T2 is varied as noted in the table. Various relaxation times and temperatures for all three experiments are 
provided. The "sim" superscript denotes a result computed directly from the simulation, and the "theo" superscript denotes a 
result calculated using the theory described in the text. 



the observed kinetic energy of the Brownian particle expressed as an effective temperature. The theory discussed thus 
far predicts that Tavg will be equal to T^s only if all the coupling coefficients are equal. Indeed, the corresponding 
temperatures listed in Table |T] are unequal indicating that a naive equilibrium treatment does not suffice. 

The effective temperature for these nonequilibrium systems can be calculated with the help of the central expression 
from Eq. (|2.29p reduced to the form of Eq. p.l2p as appropriate for the current model. The friction coefficients for 
each reservoir are an input to this equation, but they can readily be calculated in this case. They are proportional to 
the average of the reciprocal relaxation times of Eq. (|3.22p weighted by the velocity distribution f{V) of the heavy 
particle: 

2 

9klok = ^7r(i? + TkfpkUk , (3.25a) 

where 

IT = m - 41)) = / Uk{V)f{V)(fV . (3.25b) 



If indeed the heavy mass motion can be described as having an effective temperature T^s, then its velocity is simply 
the the Maxwellian distribution at the temperature Toft that is to be determined. Integration of Eq. (I3.25bp leads to 



Substitution of Eqs. (|3.25|) and (|3.26p in (|3.12p gives the self-consistent equation, 



EkV^^^Jl_ (3.27) 

The effective temperatures T'^®" calculated from this theory either by solving this self-consistent equation or by en- 
tering the inverse of the observed relaxation times into the basic Equation (j3.12p agree within the digits of accuracy 
shown in Table U and are thus included as a single value for each gedenken experiment. The good agreement between 
the effective temperatures predicted by the nonequilibrium theory and the simulated effective temperatures demon- 
strates that the coupling coefficients play a role in determining the apparent temperature of a probe particle in the 
nonequilibrium situation that it is in contact with distinct reservoirs of varying temperatures. 

The VACF relaxation time Ty can also be determined using the arguments employed thus far. It takes the value, 

2 

^ = ghok = 3 X! + TkfpkUk , (3.28) 

fe k 

which can be readily computed using the results from Eqs. (|3.25p and p.26p . The values are listed in the table as 
Ty^°°. Again, the agreement between the nonequilibrium theory and the simulation is remarkable. As a check on 
the quality of the simulations, the relaxation time te in the kinetic energy of the heavy particle was also obtained 
by simulation. The expected result form the theory of Brownian motion is that 2te will be equal to the velocity 
relaxation time as is indeed observed within nominal error bars. 



17 



4- Local Quasi- Equilibrium Limit 



In Ref. 48,, slowly relaxing systems — effectively at local quasi-equilibrium — are described on the basis of On- 
sager's fluctuation theory. Therein, the non-Markovian Fokker-Planck equation for the conditional probability density, 
P(d, t; cJg, to), was shown to be 



dP 
'dt 



d 



2J 



Xj{a, t)P 



kBTjt) dP 
m dai 



(3.29) 



where m is the particle's mass; d = {cti} is the set of fluctuating variables and do are their values at time to. The 
thermodynamic forces, Xi, are defined through the quasi-equilibrium probability density Pqo as follows: 



d 

Xi{d,t) = -TT- 



kBT{t) 



In Pqo (a, t) 



(3.30) 



The coefficients, P^ , define the connection between the stream velocities, Va-, and the conjugate thermodynamic 
forces, 



m oaj 



(3.31) 



in accordance with the rules of nonequilibrium thermodynamics. The temperature T{t) of the system is assumed to 
be time-dependent, but at a slow enough rate that the system achieves local quasi-equilibrium. This temperature is 
a function only of time and is expressed through that of the bath, TB(t), as 



T{t) ^ A{t)TB{t) 



(3.32) 



The factor A(t) is defined via the thermodynamic functions taking into account the energy exchange between the 
system and the bath and is equal to unity when the process is reversible. 

Note that the factors A{t) and Bijia, t, to) are not specified within the approach of Santamaria-Honeck et al^^ and 
other methods must be used to obtain them. Thus, the Brownian motion in a granular gas has been investigated 
using a kinetic theory to calculate these factors. ®° In this case, the behavior of a particle is described through its 
velocity, so that the stochastic variable has only one component, ai = q. The factor A has been found to be constant 
and identified through the restitution coefficient A as A = (1 + A)/2, and the single matrix element Pn = Aj{t,to), 
where 7(t, to) is the friction coefficient which depends only on time. 

The quasiequilibrium distribution is 



exp 



mq 



2fcBT(t) 



and the Fokker-Planck equation p.29p takes the form 

dP d 



dt dq 



[Bn{t,to)qP] 



d_ 



2 r 



Bii{t, to) 



kBT{t) 



It corresponds to the iLE of the general form (|3.13p 

q = -Bn{t,to)q + m 
where the stochastic force correlation function is 

Pii(t, to) 



mm) 



kBT{t)5{t~-t') 



(3.33) 



(3.34) 



(3.35) 



(3.36) 



Thus we, once again, find agreement between the current framework and earlier work. 
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IV. CONCLUSION 



In this article, we have presented a generahzed construction for the effective temperature of a tagged particle 
connected to an arbitrary number of time-dependent inhomogeneous reservoirs. It is in agreement with several limiting 
cases described earlier by various authors i^^i^^'^^i'^^'^'^i'^^ In Ref. [s^, it was shown by comparison between theory and 
simulations that the Brownian diffusion of a tagged particle (or probe) within swelling hard spheres at constant 
temperature can be surmised by an iLE. The latter is the memoryless limit of the iGLE, when the memory kernel 
^{t,t') is proportional to the (5-function, 5{t — t'). The formalism presented here is suitable for describing situations 
when the environment solvating a tagged particle is itself out of equilibrium. The particle, thus, "experiences" different 
environments with differing properties (say, temperature and viscosity) as it moves in time. 

The numerical simulation of a particle diffusing through a gas that is somehow composed of three distinct tem- 
perature particle baths discussed in Sec. IIII C 31 is instructive. Clearly, in the limit that the three baths can couple 
to each other, they will equilibrate to the same average temperature. The chosen particle will likewise attain this 
equilibrium temperature which is also equal to the effective temperature defined by Eq. (j2.29p in this limit. The 
numerical gedenken experiment provides the result for the opposite limit in which the three baths are completely 
uncoupled to each other. The perhaps surprising result to a naive observer, who incorrectly assumes that this is an 
equilibrium process, is that the chosen particle does not exhibit dynamics at the average temperature of the three 
baths. Instead, the degree of coupling between the baths and the chosen particle modulates the heat transfer that goes 
between the distinct baths through the particle. The effective temperature exhibited by the chosen particle, as given 
by Eq. ()2.29p . is precisely the mathematical form for this delicate balance. The existence of energy flows through the 
chosen particle would in the long time limit lead to the re-equilibration of all the baths. However, when the baths 
are large enough, such energy transfer is sufhciently small that it doesn't change the temperatures of the baths. Thus 
the presence of a subsystem at an effective temperature unequal to the average temperature of its nonequilibrium 
surroundings will be seen whenever the re-equilibration time of the surroundings is much longer than the time scales 
of interest in the subsystem. 

Such behavior should be seen in solution chemistry in which the solvent molecules undergo chemical reactions, 
leading to a substantial (and possibly heterogeneous) change in the solvation of the solutes. The differences between 
these neighborhoods depend strongly on the differences between the reactant and product solvent molecules. Moreover, 
the energetics of the reactions may also lead to localized energy losses or gains that would manifest themselves as 
temperature fluctuations. Although each of the solvent reservoirs would not be as simply decomposable as in the 
gedenken simulations of Sec lIII C"3l they would nevertheless interact with the solute heterogeneously leading to an 
effective temperature as given by Eq. (|2.29p . 

Another illustration lies in the nonequilibrium dynamics of colloidal microgel particles. The latter can change 
their volume in response to the temperature or pH alteration of the solutioui^ In accordance with recent unpublished 
experiments by Lyonj^ these colloids, being initially in a glassy state, can form a liquid after decreasing their size or 
crystallize when swelling. The generalization of the iGLE to incorporate temporal temperature changes thus extends 
our earlier theory^ to account for the diffusion of particles in temperature-dependent colloidal suspensions. 



This work has been partially supported by the National Science Foundation through Grant Nos. NSF 02-123320 
and 04-43564. RH is the Goizucta Foundation Junior Professor. 



In this appendix, we sketch the derivation of the solution of Eq. (|2.20p quoted in the text as Eq. (|2.23p . For 
simplicity, all the indices associated with the various baths are omitted, and consequently Eq. (|2.20p takes the form: 
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APPENDIX A 



X + v{t)x + VL^{t)x = c 



(Al) 



After the replacement, 



X = uw , 



(A2) 



the previous equation becomes 



ivu + 2uw + i'{t)wu + w[u + v{t)u + ri'^{t)u] 



cg'{t)h{t) 



(A3) 
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which greatly simphfies — cf. Eq. (jA16p — if only u satisfies the auxiliary differential equation, 

u + i^{t)u + n'^{t)u^O , (A4) 

for the initial conditions, 

m(0) 1 , u{0) = . (A5) 



The auxiliary function u in Eq. IA4[) - (|A5[) can be obtained using the steepest descent approximation. The usual 
substitution, 

u = exp(S') , (A6) 

leads to 

S^ + n^{t)^-i^{t)S-S . (A7) 

Since the right hand side (RHS) is a small quantity {v and S vanish if /i and h are constant), the zeroth-order 
approximation is determined entirely by f2: 

5*0 ±in{t) . (A8) 
Substituting Sq into the RHS of Eq. (|A7p , we obtain a differential equation for the first-order correction, S'l , 

Sl + n^t)^-iy{t)So-Sa. (A9) 
Using the definitions (|2.2ip of h' and il, the solution for S'l can be written, after some algebra, as 



S,it) = ±^^](^) Jl± ^Aln(^(t)/iW);. (AlO) 



In the case of an adiabatic change of parameters fj,{t) and h{t), the inequality 

d 





ii h 




fj, h 



< VL{t) (All) 



is satisfied. Therefore, 



and 



5i(t) « ±im - In {Kmt)) , (A12) 



Si(t) w ±i / n(t')dt' - In yjn{t)h{t) . (A13) 



The two complex solutions for u — each found by substitution of Eq. (jA6p by Eq. IA13I — can now be combined to 
ensure that the initial conditions in Eq. (jASP are satisfied, leading to 

CO&bJT{t) ( \T A\ 

where 

r(i)^i f n{t')dt'^ f^.dt', (A15) 
^ Jo Jo M(^ ) 

and the approximation is satisfied according to the inequality in Eq. (lAlip . This now leads to the desired intermediate 
differential equation in the unknown auxiliary function w, 

.. >(i) cg'{t)hit) 
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with the initial conditions, 

w{0) = a;(0) , w(0) = i;(0) , (A17) 



which follow from Eqs. (jA2p and (jA5|) . 

Eq. (|A16p can be manipulated further using the substitution, 

w = ab , (A18) 

leading to a new differential equation, 



ab - 



a 



b+ { + ) b 



This expression simplifies if only b is chosen to satisfy the auxiliary differential equation, 



^^'(^)'^(^W. (A19) 



6+ (Ki) + 2^)6 = 0, (A20) 



which gives rise to the solution, 

m 



m = TZTTT^ • (A21) 



After substitution, Eq. (jA19[) reduces to 

cg'{t)h{t) 

The solution of Eq. ((J22|l is 



u{t)q{t) . (A22) 



ait) - a(0) + T-^ / g'{s)h{s)u{s)q{s)ds . (A23) 
^(0) Jo 

The result for w can now be obtained by multiplying Eqs. (|A2ip and (|A23p . leading to 
From the initial conditions in Eq. (jA17p . it follows that x{0) = w{0) = a(0)5(0). Hence, 



The desired solution of Eq. (|2.20p quoted in the primary text as Eq (|2.23p now follows form Eqs. (jA2[) . (|A14p and 
(jA25|. 

APPENDIX B 

In this appendix, we rederive the iGLE projection shown earlier in Ref. H for the extended case of multiple 
disconnected baths that is of interest to this work. 
Substitution of Eq. ((2?23)) into Eq. (I2.19dp gives 

ft 



{u-^it'Mt'W 

ft ^fi ft' 

J2c?gi{t)h,it)u,{t) ^^.(^^,^^^^^,^^2 I g',{s)h,{s)u,{s)q{s)ds (Bib) 



dViq) ^f£i^)\_^,v^,(,(.),,) (Blc) 



dq ^ \ uji J dq 
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where we use the notation i = {i,k) as in the text, and the summation over i includes all the bath modes in all the 
reservoirs, i.e., = The second line [Eq. (jBlbp ] can be rewritten as 

X ^Y.''i9mhk{t)ui{t) J\',{s)hk{s)ui{s)q{s)ds iu,{t')!^k{t'W ^^^^ 

after a change in the order of integration, i.e., by noting that J* dt' ds = Jq ds J^dt'. The inner integral — 
appearing also in the first line (|BlaP — can be found readily with the help of the substitution dr^ = {hk{t') / ^J,k{t')) dt' 
[cf. Eq. ([223]: 

Js {ui{t')^ik{t')y cos^(wir^) cos^(wiT^) uji 

Insertion of this result in Eq. (|B2p and substitution of Ui(t) according to Eq. (|2.24p gives 

cos UJlTk(t) 



X = Y.''lg'k{t)hk{t) 



^/ tJ.k{t)hk{t) 



I coswiTfe(s) tana;irfe(t) - tanwiTfe(s) 

9k{s)hk[s) q{s) ds (B4) 

\/ fJ-kisjhkis) 



s-^ c?gkit) hk{t) f .hkis) 

2^- 7TT / 9k[s) — i-^q{s)smuJi{Tk(t)-Tk{s))ds 

y^-T^7r5fcW / gk{s)qis)d COS iViiTkit) ~ Tkis)) , 
^ uJtMt) Jo 



where g'kit) has also been expanded according to Eq. (|2.12p . Integration by parts, leads to 

^ = -E4^5feW rcosc.i(T,(i)-rfe(s)).gfe(s)g(sKs + Vf^^yg (B5) 
^ wfAife(^) Jo i V ^^i / 

_ gk{t)q{0)cosujiTkit) -Y] \^''^M 9k{t) [ cos LJi{Tk{t) - Tkis))gkis)qis)ds . 



Substitution of the results of Eqs. (jB3p and (|B5P into the differential equation (|B1[) for 5 leads to 

^^^'^^ l'\{t,t')q(t')dt'+at) 



dq 

- T^^P^Mt) [\osu;,{rk{t)~rkis))gkis)qis)ds~j^SV2iqi-),t) (B6) 
^ujffikit) Jo dq{t) 

j{t,t') = J2^!^g,{t)gk{t')cosu;,{rk{t)~rk{t')), (B7) 

^(t) = Y] —^^9k{t) (pi{0)smoJiTk + (ujiXiiO) - -<z(0) ) cosc^iTfc ) . (B8) 

Taking the derivative d{5V2) / dq with the help of the Euler-Lagrange variational principle and ignoring the higher 
order contributions from the nonstationary memory correction, as suggested by Ref. 35, one sees that this derivative 
coincides with the formal differentiation of Eq. (|2.4p in q{t), 

^ dV2 = - f dt'a{t,t')q{t') , (B9) 



dqit) 

with a(t,t') from Eq. (|2.13p . and two last terms in Eq. (|B6p cancel each other. 
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